Solution of a model of self-avoiding walks with multiple monomers per site on the 

Bethe lattice 
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We solve a model of self-avoiding walks with up to two monomers per site on the Bethe lattice. 
This model, inspired in the Domb- Joyce model, was recently proposed to describe the collapse 
transition observed in interacting polymers When immediate self-reversals are allowed (RA 

model), the solution displays a phase diagram with a polymerized phase and a non-polymerized 
phase, separated by a phase transition which is of first order for a non-vanishing statistical weight 
of doubly occupied sites. If the configurations are restricted forbidding immediate self-reversals (RF 
model) , a richer phase diagram with two distinct polymerized phases is found, displaying a tricritical 
point and a critical endpoint. 

PACS numbers: 05.40.Fb,05.70.Fh,61.41.-|-e 



I. INTRODUCTION 

Since the pioneering work by Flory fH self- and mutu- 
ally avoiding walks on lattices have been used extensively 
to model linear polymer molecules. The realization that 
the graphs relevant in the high-temperature series expan- 
sion of the n-vector model of magnetism correspond to 
self-avoiding walks in the limit n — > Q made it possible 
to apply the ideas of the renormalization group to poly- 
mer systems, and the concept of scaling has been very 
important in the later development of this field Q . 

If the polymer chains are placed in a solvent, as the 
temperature is lowered the chains may change from an 
extended to a collapsed configuration. The temperature 
where this transition happens was named the 9 tempera- 
ture [ij . This phenomenon was modelled by self-avoiding 
walks with an attractive interaction between monomers 
on first neighbor sites of the lattice (interacting SAW's - 
IS AW) , and if this problem is considered in an ensemble 
which is grand-canonical with respect to the number of 
monomers the 0-point of polymer collapse was found to 
be a tricritical point [1]. Thus, in these generalized pa- 
rameter space, the transition between a non-polymerized 
and a polymerized phase is continuous at high temper- 
atures, but becomes of first order as the temperature is 
lowered, and a tricritical point separates these two tran- 
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sition lines. The polymer-solvent system may also be 
treated in more detail, using a lattice gas model with 
two different particles (polymer monomers and solvent 
molecules) [Bj , and such models describe quite well exper- 
imental results of equilibrium polymerization in solution 
y. Actually the simpler IS AW model may be obtained 
as a limiting case of the more general diluted polymer 
model 0, Q , and the tricritical points in the two mod- 
els apparently belong to the same universality class. In 
two dimensions, extensive numerical work has been done 
to calculate the tricritical exponents of the ISAW model 
0, H, , and for the diluted polymer model the tricrit- 
ical exponents were calculated exactly [ll|, [l^. Also, it 
was found that in some cases the ISAW model may dis- 
play a third phase, besides the non-polymerized and the 
polymerized phases, which is a dense phase. This behav- 
ior seems to be restricted to two-dimensional lattices, and 
was found on the Husimi lattice [H, [3, [l^, [l^ as well as 
on the square lattice, using transfer matrix and finite size 
scaling techniques [17| . From these studies one may con- 
clude that even qualitative aspects of the thermodynamic 
behavior of interacting polymers may be determined by 
details of the model. On the square and the q = A Husimi 
lattice, for example, a dense phase exists if the interac- 
tion is assumed to be between first- neighbor bonds, but 
is not stable if interacting monomers are considered. 

Including interactions between SAW's is a major com- 
plication of the model, both in closed form approxima- 
tions and in transfer matrix calculations. In a recent 
paper [18|], Krawczyk et. al. proposed and alterna- 
tive approach which seems to be very promising, since 
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only one-site interactions are present. In this model, 
Ui — 0,1, . . . , K monomers may be placed on site i of 
the lattice, and each site contributes with a factor LUm (i) 
to the statistical weight of a configuration. Without loss 
of generalization, we may fix the weight of an empty site 
to be equal to one. This model is inspired on a gener- 
alization of a model of weighted random walks proposed 
by Domb and Joyce some time ago .19] and a it was also 
discussed in [23| for the particular choice = t^^""' of 
the weight functions. The model was studied in us- 
ing simulations for the particular case if = 3 and as a 
function of W2 and L03 , fixing uji — I. With a fixed value 
for uji only polymerized phases will appear, and the au- 
thors look for transitions between those phases. For the 
particular case of a model where immediate self-reversals 
are forbidden, simulations on a cubic lattice (model RF3) 
show a transition between extended and collapsed phases. 
As the parameters are changed, there are indications in 
the simulations that the transition changes from second 
to first order. In another case (RA2-immediate reversals 
allowed, simulations on the square lattice), no evidence 
of a phase transition is found. 



In this paper, we solve the model proposed in [T^ on a 
Bethe lattice with coordination number q. Although the 
critical exponents are classical for Bethe lattice solutions, 
they may provide better information on the thermody- 
namic behavior of models than simpler mean-field ap- 
proximations, particularly for polymer models, since first 
neighbor correlations are included [2l[. For the sake of 
simplicity, we restrict our calculation to the case K — 2, 
where collapse may already occur, but we will not as- 
sume wi = 1, so that transitions between polymerized 
and non-polymerized phases may be found. Similar to 
what happens in the studies of the ISAW model and in 
the simulations of the model with multiple monomers per 
site [l^, qualitative differences are found on the Bethe 
lattice solutions of the RA and RF models. In the more 
restrictive RF model, we find two distinct polymerized 
phases, one of them with double occupied sites only. This 
phase does not appear in the RA model. For g > 2, no 
qualitative differences in the phase diagrams are found 
as q is changed. Although a tricritical point is found in 
the RF model, as expected, in the RA model the tran- 
sition between the non-polymerized and the polymerized 
phases is always discontinuous for non vanishing values 

of UJ2- 



In section |TI] the model is defined in greater detail and 
its solution on the Bethe lattice in terms of recursion rela- 
tions is presented. The thermodynamic properties of the 
model are obtained from the recursion relations and pre- 
sented in section IIIIl Final conclusions and discussions 
are in section IIVI 



II. DEFINITION OF THE MODEL AND 
SOLUTION IN TERMS OF RECURSION 
RELATIONS 

We consider linear polymers, formed by monomers 
linked by bonds. The monomers are placed on the sites 
of a lattice and the chains are self- and mutually avoid- 
ing walks on this lattice. Up to K monomers may be 
located at the same lattice site. All monomer are sup- 
posed to be distinguishable. The polymeric chains are 
linear self- and mutually-avoiding walks, each step link- 
ing two monomers located on first-neighbor sites. Thus 
multiple steps of the walk may be placed on the same 
bond of the lattice. The chains are not allowed to form 
rings. The statistical weight of a configuration will be 
a product over site configuration weights <jj„;(z), where 
i labels the site and rii = 0,1, ... ,K is the number of 
monomers located on site i. We assume the weight of an 
empty site to be unitary (wo(i) = 1) and the weights to 
be independent of the site (u;„;(i) = The model 

may be interpreted as an interacting polymer model, so 
that uji — z — exp^fi/ksT) is the fugacity of a monomer 
and ujn = exp(— e(n)/fcBr), e„ being the interac- 
tion energy of the set of n monomers located on the 
same site. Monomers on different sites do not interact. 
Usually, if we are interested in the collapse transition, 
the monomer-monomer interactions should be attractive, 
and thus e(n) < 0. For K = 1 the usual polymeriza- 
tion model in the grand-canonical ensemble is recovered, 
which has been much studied as a model for equilib- 
rium polymerization Q . The model described here is the 
same studied by Krawczyk et al ^1^, with a slight gen- 
eralization, since they fixed 2; = 1 in their simulations. 
With the parametrization adopted by them, the poly- 
merization transition between a non-polymerized and a 
polymerized phase is not found. Following Krawczyk et 
al [iSj, we consider two versions of the model. In the 
RA (reversion allowed) variant, all walks which obey the 
self-avoidance restraint are included, whereas in the RF 
(reversal forbidden) variant, only walks without immedi- 
ate self-reversals are allowed. In other words, in the RF 
case configurations such that the walk reaches a site and 
immediately leaves this site through the same bond are 
prohibited. 

We will here discuss the solution of the model for the 
particular case K = 2 on the Bethe lattice. This is the 
smallest value for K for which the collapse may be found, 
so we have chosen it for simplicity. The Bethe lattice is 
the core of a Cayley tree, and the exact solution of sta- 
tistical models on it may be viewed as an approximation 
of the solutions of the same models on regular lattices 
with the same coordination number [22|. We consider a 
Cayley tree of coordination number q and place the end- 
points of the walks on the surface of the tree. In Fig. [1] 
a configuration which contributes to the partition func- 
tion of the RA model is shown for a tree with q — 4: and 
three generations. This walk would not be considered 
in the RF model, since two immediate self-reversals are 
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FIG. 1: Example of a configuration of the Cayley tree with 
q = 4 and 3 generations. This configuration is allowed in the 
RA model, but not in the RF model. Polymer bonds (steps 
of the walks) are represented by full lines, while the lattice 
bonds are dashed lines. 




FIG. 2: Contributions to the recursion relation for the partial 
partition function gi (RF model). 



located at the new site. The recursion relation for is: 

(la) 



I q-1 

9o,i = % ' 



9o,2 



2 



X 1 t^i5o ^9l, 



5o,5 



The terms of the recursion relation for g[ are: 
g'l.i ^ (9-1) X l^i5o"^5i, 



5l,2 



3 



For g'2 we find two terms: 
'q-1 



52,2 = (9 - 1) X 2 W2 5o ^92 



(lb) 
(Ic) 



9oA = ^('^ -s^) x4c.2ffo'''5i52, (Id) 



(le) 

(2a) 
(2b) 



5i,3 = 2 ^ ) X 4^250 5152. (2c) 



(3a) 
(3b) 



present. The statistical weight of the depicted configura- 
tion is oj\lo2- 

The solution of the model on the Bethe lattice is ob- 
tained by defining partial partition functions on rooted 
sub-trees, and realizing that a subtree with n+\ gener- 
ations may be obtained through the connection of g — 1 
n-generations subtrees to a new site and bond [2ll | . This 
operation leads to recursion relations for the partial par- 
tition functions. We start with the more restrictive RF 
model. One should be careful with the multiplicity of 
each contribution to the recursion relations, recalling 
that even monomers placed on the same site of the lattice 
should be considered distinguishable. The partial parti- 
tion functions for the model are (70 j 5ij and 32, where the 
subscript denotes the number of polymer bonds on the 
root bond of the subtree. In Fig. [5] the contributions to 
the recursion relation for gi are shown graphically, in the 
order they appear in the equations [2l 

Below the recursion relations are presented. In gen- 
eral, we have g[ — J2j9ij^ and, whenever appropriate, 
the terms in the sums begin with a product of two nu- 
merical factors, the first of which is the multiplicity of 
the configuration of the incoming bonds and the second 
is the multiplicity of the connections with the monomers 



We may now define ratios of the partial partition func- 
tions Ri — 51 /go and i?2 — 52/50- Any expected value in 
the center of the tree may be expressed as a function of 
these ratios, so that if we are interested in the behavior 
of the model in the center of the tree and in the ther- 
modynamic limit, we should find the fixed point of the 
recursion relations for the ratios, which are: 



where 

D = 1 



(g - l)cji + 6( ^ ^ "'^ ) CJ2-R? 



UJ2 i?2 



D ' 



9 - 1 \ n2 



Ri + {q- 1) i?2 



(4a) 



9' 



^^\l02RIR2+ ^{^ ^^']0J2Rl (5) 



For the RA model, where immediate self-reversals are 
allowed, three other root configurations for subtrees are 
possible. The additional partial partition functions are: 
f/3 (two polymer bonds at the root, connected to the same 
monomer above) , 54 (same as for , but with more than 
one monomer in the closed path), and 35 (three polymer 
bonds at the root, two of them connected to the same 
monomer above) . The recursion relations for this gener- 
alized model are: 



to g[ are: 



5o,i 
5o,2 



50,3 



5o,4 



5o,5 



fQ,6 



5o,7 



50,9 



fo.ii 



9-1 



^ ^ 9-3 2 

2 ) X 1 ^1 5o 5i , 

^^^\ ^ r 9-5 4 
^ I X 6 tJ2 5o 9i, 



3 ( 2 ) X 6 W2ffo ^ 9l92: 



? ^ 1 \ ^ r 9-3 2 
2 I X 6tJ2 5o 52, 



(g - 1) X 1 ujigl ^32, 

X 2 cjaffo^'^Si 53, 



3 



g-1 
3 

2 
2 
2 



X 4 cj2 5o "51 54, 



X 2 W2 5o ^ 52 53 



X 4 cj2 5o ^5254, 



X 2 W2 5o ^51 55 



(6a) 

(6b) 

(6c) 

(6d) 

(6e) 
(6f) 

(6g) 

(6h) 
(6i) 

(6j) 
(6k) 



5m 



5l,2 



5'l,3 



5l,4 



5U 



(9 - 1) X 1 uigl ^gi, 



3 



r- 9 — 4 3 

X 6 W2 5o 5i 



2(^2^) X 6 tJ25o ^5152, 



2( 2 M X 2^250 '5153, 
2(''^ 2 x4^2 5o"'5i54, 



5^,6 = (9 - 1) X 2cj2 5g ^55, 



(7a) 

(7b) 

(7c) 

(7d) 

(7e) 
(7f) 



The contributions to 172 identical to the ones in the RF 
model, shown in equations [H For 33 we find: 



53,1 



53,2 



53,3 



9-1 
^i5o > 



\^)x2u:,gl-^gl 



(9 - 1) X 2 W2 5o ^52 



The contributions to 174 are: 

54.1 = (9- 1) X 1 ^2 3^"^ 33, 

54.2 = (g - 1) X 2 W2 ffo"^ 34- 

Finally, a single contribution to g'^ was found: 



(8a) 

(8b) 
(8c) 



(9a) 
(9b) 



55,1 = (9 - 1) X 2tj2 5o ^51 



(10) 



Notice that, although the first five terms are similar to 
the ones in the RF model, the multiplicities related to 
the connections of the two additional monomers may be 
larger in the less restrictive RA model. The contributions 



In the recursion relations for the RA model, we notice 
that the partial partition functions 173 and g^ always ap- 
pear in the combination g^-\-2 g^. We therefore define the 
following ratios for this model: i?i = 51 /go, R2 — 52/50, 
R3 = (53 + 2g4)/5o, and R4 = 55/50- The recursion 
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relations for the ratios are: 
R[ = 



12( 2 M LU2R1R2 +4r 2 ^ ) ^2RiRz 



R'2 = 2 



2{q-l)uj2RA 

g-1 



1 

15 



Rf + {q~ 1) i?2 



W2 



2(g-l)cj2 (i?2 + i?3) 

, 2 (g - 1) Ri 
R4 = ^ , 



1 

D 



(11a) 
(lib) 

(11c) 
(lid) 



where 
D = 1 



\^y,Rl+ 6(\')u;2Rt 



{q-l)uJiR2 + 6(^'^ 3 cc;2i??i?3 + 
4('^2 LO2R2R3 +4r 2 j i?i -^4. (12) 



The partition function of the model on the Cayley tree 
may be obtained if we consider the operation of attaching 
q subtrees to the central site of the lattice. The result for 
the RF model will be: 



12 ( y gl' ff? 52 + 4 r j C02 dl (13) 



For the RA model, a similar calculation leads to the par- 
tition function: 

Y^^^^ = gl + (^) gr' gl + 6 (^) .2 gV + 
18 c^2 5o"^ 3? 52 +6 UJ2 gT'^ gl + 
9 '^150^^52 + 6^^^ '^2gT^ gl{.g-i + 254) + 

^ ^2 5o"^ 52 (.93 + 2.94) + 



9-2 



2)^^255 5155 



(14) 



In the thermodynamic limit, the solution of a model on 
the Cayley tree usually shows a behavior which is quite 
different from the one expected on regular lattices, since 
the number of surface sites represents a non-zero frac- 
tion of the total number of sites, even in the thermody- 
namic limit. We therefore study mean values calculated 
at the central site of the tree. The behavior of a model in 
the thermodynamic limit and in the central region of the 
Cayley tree has been named the Bethe lattice solution of 
this model [22|. Using the partition functions above, we 
then proceed calculating the densities at the central site 
of the tree. The density of monomers is given by: 



P = 



P + 2Q 

1 + p + g' 



where for the RF model we have: 



R\i 



(15) 



(16a) 



Q 



{RF) 



UJ2 



Q\^.\uj2R\ + 12( l) RIR2 + 



Ri 



For the RA model, we have the results: 

p{RA) 



Q 



{RA) _ 



UJ2 



l]Rl + qR2 



(16b) 



(17a) 



181 ^ji?? i?2 + 



6(3ji??i?3- 



4(2)^2i?3 



4( 2 I i?i ^4 



(17b) 



It is supposed in the expressions above for the density 
of monomers at the central site that the ratios Ri have 
their fixed point values, so that the thermodynamic limit 
results are obtained. Since P and Q are non-negative, 
the density of monomers will be in the interval [0,2], as 
expected. The probability that the central site is occu- 
pied by a single monomer will be pi = -P/(l -I- P -I- Q), 
and the probability to find two monomers at the central 
site is p2 = Q/{l + P + Q). 



III. THERMODYNAMIC PROPERTIES OF THE 
MODEL 

As discussed above, the thermodynamic properties of 
the model on the Bethe lattice are determined by the 
fixed point values of the ratios, defined by the recursion 
relations, for fixed values of the parameters toi and uj2- 
For the RF model, an inspection of the recursion rela- 
tions m shows that a fixed point with Ri = R2 = may 



exist. Equations [TSl and [TBI show that in such a phase the 
density of monomers vanishes, so that we may recognize 
it as the non-polymerized (NP) phase. The next step 
is to identify the region of the parameter space (0^1,072) 
where this phase is stable. We thus calculate the jacobian 
in the NP phase, whose elements are: 



J, 



dR' 



(18) 



and the NP phase will be stable in the region of the pa- 
rameter phase where the absolute values of both eigen- 
values of the jacobian are smaller than 1. The jaco- 
bian matrix is diagonal and the result of this calculation 
shows that the NP phase is stable when the conditions 
uji < l/{q — 1) and W2 < l/[2(g — 1)] are both satisfied. 

Another fixed point of the recursion relations for the 
RF model is such that i?i = and i?2 7^ 0. In this phase 
pi vanishes but p2 not, so it is a polymerized phase in 
which only double occupied (DO) sites are present. The 
fixed point value for Ri is given by: 



R 



DO 



2((?- l)c^2 - 1 



2(9-l)(g-2)c^2' 
and from [16] we have P = and 

Q = 4[2((7-l)a>2-l], 
which leads to the density of doubly occupied sites: 
i(g- l)t^2 - 4 



(19) 



8(q- l)c^2 



(20) 



The jacobian matrix is also diagonal in the DO phase, 
the non-zero elements are : 



Ji,i = 



J^a — 



2UJ2 



l2{q~2) [2(g-l)u;2 - 1] 



{q - 1)^2 



1 



{q — 1) UJ2 



,(21a) 
(21b) 



The DO phase is stable if the absolute value of both 
eigenvalues of the jacobian is less or equal to one, which 
leads to the conditions: 



1 



<LlJ2< 



2{q-2) 



2(9-1) - - (g-l)(4(7-9)' 
and for a value of UJ2 in the range above, we have: 



(22) 



^1(^2) < 2ti)2 



12 [2(g-l)c^2 - 1] {q-2)uj2 



(9-1) 



It may be observed that this phase has a non-vanishing 
region of stability for any value of the lattice coordination 
number. For large values of q we may expand the limits 
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FIG. 3: Phase diagram for the RF model on a Bethe lat- 
tice with g = 4. Full lines are continuous transitions and on 
the broken lines two phases coexist. NP denotes the non- 
polymerized phase (pi = p2 = 0), DO the double occupancy 
phase (pi = 0, p2 7^ 0, and P the regular polymerized phase 

(Pl,P2 /O). 



of the stability interval in powers of 1/q, finding that: 



1 



2(9-1) 



2q \ q 



2(9-2) 



(9-l)(49-9) 



1 

2^ 



1 



We conclude that the NP-DO phase transition is contin- 
uous, ati:j2 = l/[2(q-l)]. For W2 > l/[2(g-l)] the other 
stability limit of the DO phase is given by expression [23l 
A third fixed point of the recursion relations [31 which 
we may associate to the regular polymerized phase (P), 
corresponds to both ratios at non-zero values. The an- 
alytic study of this fixed point is harder, but with the 
help of a formal algebra software we may conclude that 
the stability limit of the DO phase is never coincident 
with the stability limit of the P phase, and thus the 
DO-P transition is always of first order. The stabil- 
ity limits of the NP and the P phase are the same for 
u!2 < l/[2(9 — 1)^], henceforth the transition between 
these phases is continuous in this region. For larger val- 
ues of UJ2, this transition becomes discontinuous. A tri- 
critical point is found at: 



TC 



TC 



1 



9-1' 
1 



2(9-1) 



(24a) 
(24b) 



The first-order transition lines may be found through 
a Maxwell construction [l3| , but this procedure may be 
difficult to converge numerically if more than two phases 
are present. For simplicity, we thus adopted an alter- 
native procedure (l4j] . iterating the recursion relations 
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0.20 



FIG. 4: Densities for the RF model with g = 4 as functions 
of u>2, for uJi = 0.09. The full line is the density of double 
occupied sites (p2) and the dashed line corresponds to the 
density of single occupied sites (pi). 



starting with the physical values for the partial partition 
functions for a subtree of generation "zero" , which are 

^o"^ = 1, g^i^ = wi, and g^^ — lu2- With these starting 
values, the recursion relations do converge to the sta- 
ble thermodynamic phase and therefore we may obtain 
the first-order transition lines. In Fig. [3] we show the 
phase diagram of the RF model for g = 4. The continu- 
ous lines are second order transitions and at the broken 
lines two distinct phases coexist. The continuous NP- 
P transition lines ends at a tricritical point located at 
= - 1) = 1/3 and UJ2 = l/2(g - 1)^ = 1/18. 
The second order NP-DO transition ends at a critical 
endpoint at lji = 0.1190365 and u!2 = 1/6 (wi was de- 
termined numerically). It is worth mentioning that there 
is no discontinuity at the point where the P-DO coexis- 
tence line touches the wi = axis, denoted by a square 
in the phase diagram. Also, the fact that the DO phase 
becomes unstable as uji is increased for fixed W2 may be 
understood noting that although it is energetically more 
favorable to have doubly occupied sites at high values of 
^2, a non-zero density of simply occupied sites increases 
the entropy. 

In Fig. m the densities of single- and double occupied 
sites are shown as functions of W2 for a fixed value uji = 
0.09, again for a lattice with q = 4. For low values of 
u!2 the NP phase is stable. At W2 = 1/2(9 ^ 1) = 1/6, a 
continuous transition to the DO phase occurs, followed 
by a discontinuous transition to the P phase at uj2 ~ 
0.17503. 

Turning our attention to the more general RA model, 
we notice from the recursion relations [Tl] that the non- 
polymerized phase has now the fixed-point values i?i = 
R2 — R4 — and 



These values for the ratios actually correspond to p = 0, 
as may be seen using equation [15] for the density of 
monomers. The stability limit of the NP phase is ob- 
tained requiring that the largest eigenvalue of the jaco- 
bian, calculated at the fixed point values of the ratios 
above, has an absolute value equal to 1. This leads to 
the following condition: 



UJ 



{NP) 



[l-2(g-l)c^2]'[l + 2(g-l)c^2] 
l-2u;2 



(26) 



We may search for a DO fixed point, which should 
have the values Ri = R4 — and i?2,^3 7^ 0. How- 
ever, writing the equation for the fixed point value of 
the ratio R2/R3, we reach the unphysical conclusion 
i?2 = —u>i/[2{q — l)a-'2]7 indicating that such a phase 
does not appear in the RA model for wi ^ 0, since ra- 
tios of partial partition functions should be non-negative. 
We may have some physical understanding of the non- 
existence of a region where the DO phase is stable in 
the RA model if we notice that the DO phase found in 
the RF model is composed by linear chains with double 
polymer bonds between first neighbor sites occupied by 
two monomers, which are constrained to have endpoints 
at the surface of the tree. We may notice that in the RF 
model, for wi = 0, we have a continuous transition from 
the NP phase to the DO phase at uj2 = l/[2(g - 1)]. If 
we take into account that at each lattice site there are 
two ways to connect the bonds to the monomers on this 
site for a linear DO polymer, this corresponds exactly to 
the well known result that the polymerization transition 
for linear chains on a Bethe lattice occurs when the sta- 
tistical weight of a monomer is equal to l/(g — 1) [2l| . 
as may be also seen in our results for the RF model in 
the particular case tt'2 = 0. If we allow immediate self- 
reversals, these chains may end at any step inside the 
tree. Therefore, this corresponds to allow for endpoints 
in the problem of polymerization of linear chains, and 
it is known that this effect eliminates the polymerization 
transition. It is possible to map the equilibrium polymer- 
ization model into the n-vector model of magnetism in 
the formal limit n — > 0, and if this is done, the statistical 
weight of monomer at the end of chains is proportional to 
the magnetic field in the n- vector model 23]. The con- 
tinuous ferromagnetic transition in the n-vector model 
exists only when no magnetic field is applied. Accord- 
ingly, no polymerization transition is found for nonzero 
statistical weight of monomers at endpoints of chains. 

The phase diagram for the RA model is shown in Fig. 
[5l The coexistence line between the NP and P phases 
was again found iterating the recursion relations 1 1 1 1 with 



initial values gg"' = 1, gj' 

.gf ^ = 0, and g^°'' ^ uj2. In the RA model, the NP-P 
transition was found to be continuous only at vanishing 
values of uj2, so that the phase diagram is qualitatively 
different from the one found for the RF model. At {uii = 
l/(q — l),uj2 ~ 0) we have, therefore, a critical point and 
no tricritical point is present in the phase diagram. 
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FIG. 5: Phase diagram for the RA model on a Bethe lattice 
with g = 4. Only the non-polymerized (NP and the regular 
polymerized (P) phases appear. For non-zero iJ2, the phase 
transition is always discontinuous. The stability limit of the 
NP phase is shown as a dash-dotted line. 



IV. CONCLUSION 

A model of self- and mutually avoiding chains which al- 
lows for up to X = 2 monomers per lattice site was solved 
on the Bethe lattice with general coordination number q. 
Endpoints of the chains are allowed only at the surface of 
the tree. Following [l^, besides the general model (RA), 
a restricted version (RF) where immediate self-reversals 
of the chains are forbidden, was also considered. Al- 
though in some respects this model resembles models for 
branched polymers which were already studied on the 
Bethe lattice f23|, the multiplicities of the contributions 
to the partition functions are different, since each config- 
uration has to be a self avoiding walk passing through all 
monomers placed on the lattice. The statistical weight of 
a configuration is lu^'^ '-^2^, where Ni {N2) is the number 
of lattice sites with 1 (2) monomers. For low values of 
uji and UJ2, a non-polymerized phase is stable for both 
models. The phase diagram for the RF model displays 
two distinct polymerized phases, one of which (DO) with 
double occupied sites only. In this model, the transition 
between the NP phase and the regular polymerized (P) 
phase may be of first or second order, and therefore a 
tricritical point is found, which could be associated with 
the 0-pomt. Also, the second-order transition between 
the NP and the DO phases ends at a critical endpoint. 
In the solution of the more general RA model only one 



polymerized phase is found, and the transitions between 
it and the NP phase is always discontinuous for nonzero 
UJ2- Thus, no tricritical point is found in this case. In 
both models no dense phase was found, similar the one 
found for the ISAW model with interacting bonds on the 
g = 4 Husimi lattice [HHIIll. 

One conclusion of this study is that both versions of 
the model have qualitatively different phase diagrams 
on the Bethe lattice. It is worth mentioning that sim- 
ulations of the model for K — 3 also presented results 
which were strongly dependent upon the restrictions of 
the model [ll], although no direct comparison can be 
done between our results and the simulations, since they 
correspond to different values of the maximum number 
of monomers per lattice site K. Although at least for the 
RF model a tricritical point is found which may be as- 
sociated to the collapse transition of the chains, if we 
consider the location of this point we are lead to an 
inconsistency. We may adopt a physically reasonable 
parametrization of the two statistical weights, stating 
that uji — z — exp(/?/i) is the fugacity of one monomer 
and UJ2 — z^LO, where w is the Boltzmann factor associ- 
ated with the interaction of the two monomers placed on 
the same site (interactions between monomers at differ- 
ent sites are not considered in the model). The tricritical 
point, whose location is given by the equations 1241 corre- 
sponds to a fugacity z — l/{q — \) and a Boltzmann fac- 
tor Lj — 1/2. Since lj < \, the collapse transition would 
occur for repulsive interactions between monomers occu- 
pying the same site. For the ISAW model on the Bethe 
lattice, the tricritical point is located at z — l/(q — 1) 
and w = {q — l)/{q — 2) > 1 Q, in the attractive re- 
gion. It should be mentioned, however, that the solution 
of the ISAW model (interaction between monomers) on 
the Husimi tree with q = A leads to a tricritical point at 
Lo « 1.54 [l3|, which is somewhat higher than the Bethe 
lattice result. It may be interesting to see if the solution 
of the RF model on the Husimi tree (which is supposed 
to be closer to the solution on regular lattices), displays 
a tricritical point in the physically expected attractive 
region. We are presently working on this problem. 
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